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Abstract 
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I. INTRODUCTION 



The phenomenon of entanglement is one of the most fundamental and counterintuitive 
features of quantum mechanics. Its fundamental role was emphasized by the formulation of 
the EPR paradox [1], despite the original purpose of the latter to question physical reality 
of the wave function. The counterintuitive nature of entanglement is a hallmark of quantum 
mechanics, and its properties reveal deep distinctions between quantum and classical objects. 

Prom the mathematical point of view entanglement is a consequence of the superposition 
principle and the tensor product postulate in quantum mechanics. Specifically, the principle 
and postulate imply that a state vector of a system consisting of several subsystems is a 
linear combination of tensor products of state vectors of the subsystems. Mathematical and 
physical properties of states interrelate: a state is disentangled if it can be transformed into 
a factorizable state; any other state is entangled. Equivalently, a state is disentangled if and 
only if each subsystem is in a definite state. 

Despite the simplicity of the above qualitative features of entanglement the complete list 
of its quantitative characteristics is unknown. For example, it might appear that the smallest 
number of linearly independent factorizable terms representing a state is an appropriate 
characteristic of its entanglement. This is true for two subsystems, in which case this single 
quantity completely classifies all entangled states. For more than two subsystems, however, 
this quantity does not characterize entanglement since it depends on a choice of bases. 
To choose appropriate entanglement characteristics for the general case we need to study 
invariant properties of states of composite systems; these are the key properties shaping the 
following discussion. 

For states of composite systems entanglement quantifies ways in which states of sub- 
systems contribute to linear combinations of tensor products. The larger the numbers of 
contributing states of subsystems, the greater the variety of arrangements of terms in linear 
combinations. However, some of the arrangements should be considered as dependent since 
they are related by transformations of bases. Such related states form equivalence classes, 
finding the complete set of which is the goal of entanglement classification. 

To classify entangled states one usually employs entanglement invariants, which are cer- 
tain invariant quantities associated with the states. The nature of the problem requires that 
the invariants do not change under all transformations that can be reduced to changes of 
bases. Consequently, the invariants take the same values for all states within each equiva- 
lence class, and the standard method of finding them uses the classical theory of invariants 
[2] . Variants of this method are used in most known cases of partial or complete entangle- 
ment classification; see, for example, [3-20]. 

Developing the ideas outlined above, we have introduced in [21] a new entanglement 
classification method based on algebraic properties of tensor products of linear maps. In 
this paper we generalize and expand both the method and its applications. We first introduce 
various equivalence relations and corresponding equivalence classes on linear spaces of states. 
We then show how these classes lead to various linear subspaces and their invariants, which 
are the central objects in our method of algebraic classification of entangled states. During 
the development of our method we turn repeatedly to the example of three qubits to illustrate 
the procedure. Finally, we proceed with numerous more complicated but physical relevant 
examples demonstrating the use of the method in classifying the entanglement of many 
systems unsolved until now. 
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II. PRELIMINARIES 



Tensor product space 

We begin by introducing the main components of our construction. Let S he a quantum 

system that consists of n subsystems {5'j}jg/, where I = {l,...,n}. For each i E I, let 
a finite-dimensional vector space Vi over a field F be the state space of Si. Extension to 
infinite-dimensional spaces is nontrivial and is not considered here. We choose F = M for 
the simplicity of presentation; the case F = C needs only minimal modifications. 

Our first task is to define the space V, the state space of S. The tensor product postulate 
in quantum mechanics says that K is a subspace of the tensor product space ^i^iVi. A 
specific choice of subspace V depends on the nature of S. For identical subsystems, for 
example, the permutation symmetry acting on the subsystems determines V. In particular, 
for bosonic or fermionic subsystems V is, respectively, the symmetric or antisymmetric part 
of the product V^. Also, if there is an equivalence relation among elements of V (as, for 
example, for linearly dependent vectors in quantum mechanics), then V is the appropriate 
quotient set. Modifications due to these and similar properties can be easily included into 
the following development, which assumes the simplest case where V — ^ie/V^. 

Transformation group 

We aim to study properties of the system S related to its composition in terms of the 
subsystems {Si}i^i; these are equivalent to properties of V related to its composition in terms 
of {Vi}ig/. The latter manifest themselves in their transformations under an appropriate 
group. Note that the tensor structure of V implies that the transformation group relevant 
for studying properties of V is not the general linear group of V, GL{V), but rather its 
subgroup Xi^jGL{Vi). Accordingly, for each i E I we choose a subgroup Gi of GL(Vi) and 
define the corresponding subgroup G = Xi^iGi of GL{V). As a result, the group G is the 
transformation group for V, and it determines properties of V related to its composition 
in terms of {Vi}i(zj. Particular cases (where only certain subsets of V and subgroups of G 
matter) are of interest as well and can be treated similarly to the general case of \^ = ^iei^i 
and G = Xig/GL(Vi) considered here. 

Equivalence classes 

The group G induces the equivalence relation ~y on V, which is given by v' ~y v for 
each v,v' eV if and only if there exists g E G such that v' — gv. The equivalence relation 
defines the equivalence class of v under ~y, 

G{v) ^{v' eV:v' ~y v}. 

Since all elements of the class C{v) are equivalent, any one of its elements determines the 
whole class. It is thus convenient to replace C{v) with its arbitrary single element v E C{v), 
which we call a representative element of the class. (For each specific class G{v) the choice 
of V based on various symmetry considerations generally leads to simplifications.) Repeating 
this procedure for each v eV, we partition V into the set of equivalence classes 

C^U,ev{C{v)} 
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such that each vector in V belongs to one and only one class. Finally, replacing each class 
in C by its representative element, we arrive at the set 

V^{ve C{v) : C{v) e C}, 
which can also be written as the quotient set V — V/ 

Properties of entangled states 

Understanding the structure of the quotient set V is our ultimate goal. We begin with a 
general property of V, its partition into three characteristic subsets of vectors: (1) the zero 
vector, (2) decomposable vectors, (3) nondecomposable vectors. By definition, a decompos- 
able vector f e y is a vector that can be written in the factorizable form v — (^lieiVi, where 
Vi & Vi is a nonzero vector for each i G /. A nondecomposable vector is a vector which 
is neither zero nor decomposable. We will derive the general form of a nondecomposable 
vector after we establish its invariant characteristics. 

The above partition is physically significant because it is in a one-to-one correspondence 
with the partition of quantum states into three types: (1) the vacuum state, (2) disentangled 
states, (3) entangled states. The zero vector (the vacuum state) and decomposable vectors 
(disentangled states) are the simplest elements of V; although they comprise only a small 
part of V, they span all of it. By contrast, nondecomposable vectors (entangled states) are 
more complex and difficult to categorize. The difficulty is combinatorial because decom- 
posable vectors from V that enter the linear combination representing a nondecomposable 
vector differ by ways in which linearly independent vectors from {Vi}i^i enter the expression. 
Finding all such possibilities of nonequivalent combinations (which is the same as finding 
the quotient set V) is the problem of entanglement classification. 

Another general property of V concerns the number of its elements. Although the set V 
is not a vector space, we use the notation dim V for the number of unconstrained elements 
of F that a general element of V depends on. Using a similar notation for dim G, we find 

dim V > dim V — dim G. 

The inequality sign appears here because, in general, the system of linear equations for 
g & G that follows from the equivalence condition v' — gv is not linearly independent. We 
have two distinct cases here: (1) if dim V — dim G < 0, the above inequality does not tell us 
if there are any unconstrained elements of F that a general element of V depends on; (2) 
if dim^ — dimG > 0, there are at least diml^ — dimG such elements of F. Consequently, 

V is an infinite set in the second case. Asymptotically for large n, dim V is exponential in 
n and dim G is at most quadratic in n. It follows that n does not need to be very large for 
the set V to be infinite; in other words, V is typically infinite. 

Example of three qubits 

To illustrate the concepts introduced above, we consider a particular example of three 
qubits, in which case there are three 2-dimensional spaces Vi,V2, V3 and their tensor product 

V — Vi<S>V2<S>V3. We choose arbitrary bases {eij}i<j<2 for each Vi and expand an arbitrary 
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element v &V in terms of its coordinates {vj^j^j^}, 

2 2 2 



jl=l j2=l j3 = l 



The transformation group G acts on v according to v ^ 919293^, where Qi e GL{Vi) and the 
coordinates of 919293V are 

222 

i9l9293v)ji,j2,j3 ^ Z Z ^fci ,fc2 ,fc3 (gl ) fci , ji (ff2 ) fca , j2 (ga) fc3 ,j3 ■ 

kl = l fe2 = l ^3=1 

The simplest example of a decomposable vector in V (a disentangled state in S) is v — 
® 62,1 ® 63^1- This is a state of the system 5" in which each of its subsystems Si is in a 
definite state Cij. 

There are three types of nondecomposable vectors in V (entangled states in S). For 
the first type, vectors are nondecomposable for the tensor product of two spaces but de- 
composable for the tensor product of the three spaces. Choosing Vi and V2 as two such 

spaces, we have the state v = (ei.i (8> 62,1 + ei_2 ® 62,2) ® 63,1 in which the subsystems 5*1 and 
5*2 are not in definite states, while the subsystem S3 is in a definite state. The other two 
states of this type are obtained by permutation of the subsystems. For the second and third 
type (which in the hterature are called respectively the W and the GHZ classes), states 
are nondecomposable for the tensor product of the three spaces. The standard forms of 
their representative states arc v = ® C2.1 ® ^3,2 + ^1,1 ® 62^2 ® ^3,1 + ^1,2 ® 62.1 CB) 63^1 and 
V = Ci^i 62,1 ® 63^1 + ei^2 ® 62,2 ® 63^2, respectively. For these states, no subsystem is in a 
definite state. 

Applying all elements of the group G to a representative vector v in any of the above three 

types of elements of V, we obtain the equivalence class C{v), which leads to 7 equivalence 
classes (counting permutations and including the zero vector which is in its own equivalence 
class). This is a well-known result (which we also proved by using our method in [21]) that 
these 7 classes constitute the complete entanglement classification of three qubits. 



Invariants 

The problem of finding V can be solved by direct or indirect methods. In a direct 
method, one uses the definition of V to derive the general form of representative elements 
of equivalence classes. Although there are no restrictions to such methods in theory, they 
are usually inefficient in practice because of the need to solve complicated equations. By 
contrast, in an indirect method, one seeks quantities characterizing elements of V which 
are invariant under G. Equivalence classes arc obtained by finding allowed values of these 
invariants. Indirect methods are usually efficient if all invariants are known. 

Continuing with indirect methods, let a{v) e F be an invariant of v induced by the group 
G. This is a quantity that satisfies a{v') — a{v) for each v eV,v' e G{v), which implies that 
invariants depend only on classes. Let A{v) be a complete set of algebraically independent 
invariants of v, so that v' ~v' v if and only if A{v') = A{v), for each v, v' G V. The standard 
method of finding A{v) is to use the classical theory of invariants and covariants; for a 
modern introduction, see, for example, [2]. Almost all known cases of partial or complete 
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entanglement classifications use this method to a certain extent; see, for example, [3-20]. 
The rapid increase of with n is the main reason why only the simplest cases of 

entanglement classification have been fully carried out. 

Let us now consider a typical case of infinite V. We find that the set of all possible values 
of the invariants, U^,gy{yl(f )}, is infinite. The resulting information about V in terms of 
its elements and invariants is both overwhelming in its detail and impractical in its use. As 
a key part of our method, we reduce the amount of information by grouping equivalence 
classes into a finite number of sets. The grouping is determined by certain equivalence 
relation between classes in each set, a natural choice for which is defined as follows. 

Equivalence of invariants 

We first introduce the rescaling equivalence of invariants. We note that since linearly 
dependent vectors in quantum mechanics correspond to the same physical state, we require 
fv e C{v) for each G V, / G F, / 7^ 0. It follows that algebraic invariants are homogeneous 
polynomials; consequently, zero is the most important value of each invariant. This suggests 
we extend the above rescaling equivalence of states to the rescaling equivalence of invariants. 
Specifically, we define the equivalence relation ~^ on the field F by setting a' r^p a for each 
a,a' E F if and only if there exists / G F, / 7^ such that a' = fa. (For F = R or F = C, 
this simply means that any two nonzero elements are equivalent.) It is easy to generalize 
this equivalence to ordered sets over F, so that for each pair of such sets {a'f,)k^K and {a^jk^K 
we define (afc)^^/^ ~f {o-klkeK if and only if a'^ ~^ for each k G K. 

Having established equivalence for invariants, wc transfer it to vectors. Namely, we define 
the equivalence relation ~y on the set V by setting v' ~y v if and only if A{v') ~f ^("w), 
for each v' G V . Since invariants depend only on classes, v' ~y v implies v' ~y v. The 
relation ~y defines the quantities C"(v), v', C", V' in the same manner as the relation ~y 
defines the quantities C(f), "C, C, V . Clearly, C is a partition of C. 

The sets C and V' are the main objects of our study. We call the problem of finding 
them the restricted entanglement classification problem to emphasize that we seek only sets 
of classes, not the classes themselves. One way to solve the problem is to use the set of 
invariants A{v) from the standard classification method. This approach requires studying 
conditions under which elements of A{v) are zero. If A{v) is known, this method gives the 
solution; however, we prefer a simpler approach that uses new algebraic invariants N{v) 
instead of A{v). The advantage of our approach is that each element of N{v) describes 
certain algebraic properties of v and takes a vahic from only a finite set of integers. The 
construction of N{v) uses only basic linear algebra [22] and proceeds as follows. 

III. METHOD 
Outline 

The set of invariants N{v) is uniquely determined by the following conditions. First, 
N{v) depends only on the equivalence class C'{v) to which v belongs. As a result, both 
C'{v) and N{v) are invariant under the action of the transformation group G. Second, the 
rescaling property of A{v) implies that N{v) depends only on properties of linear subspaces 
of V; let L{v) be the set of such subspaces. Third, L{v) depends linearly on v. Fourth, L{v) 
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describes properties of v associated with all partitions of the system 5* into subsystems build 
from {S'jjig/. Such partitions result from all choices of writing V as the tensor product of 
spaces built from {Vi}i^i. 

Maps 

The above conditions require that L{v) is defined in terms of linear maps. We find these 
as follows. We first partition the system S into subsystems T and T', so that S = T U T' . 
Let W and W be the state spaces for T and T', respectively, so that V = W ® W. Our 
main tool for constructing L{v) is a linear map 

f{v):W^W', f{v){w) = v(^w*, 

where w* e W* is the dual of w eW. 

According to a standard result in linear algebra, all information about a linear map is 
contained in two fundamental spaces associated with it: its kernel and image, 

ker/(v) = {weW: f{v){w) = 0}CW, 
imf{v) = {w' e W: w' = f{v){w), w e W} C W' . 

Associated with the map f{v) is the transpose map 

f'iy) -.W ^W, f{y) {w') =v®w'*. 

The matrices of f{v) and f'{v) are the transposes of each other. 

Introducing inner products in W and W , we can relate the kernels and images of f{v) 
and f'{v) through orthogonal compliments, 

imf{v) = (kerr(.;))^, imfiv) = {kci f{v))^. 

(The orthogonal complement Y-^ of a subspace Y of an inner product space X is the set of all 
vectors in X that are orthogonal to every vector in Y, Y^ = G X : {x,y) = 0, \/y e Y}-) 

Thus, if both maps are used to construct L{v), then it suffices to consider only their 
kernels, for example. We adopt this choice. Furthermore, since f'{v) is obtained from f{v) 
by interchanging W and W, both maps are included by considering only f{v) for both 
V (^W andV = W' ^W. 

For specific computations we need expressions for the above quantities in terms of coor- 
dinates. We introduce these by choosing arbitrary bases {ei}i<i<dimw ^-nd {e^}i<i<dimTy for 
the spaces W and W and representing a vector v e V in terms of its coordinates {vij}, 

dim W dim W 

i=l i=l 

We find 

dimWdimW 
i=l j=l 

dim W 

ker/(v) = |w e VT: ^ VijWi = 0, j e {1, dim Vr'jj. 

i=l 

The kernel of a map is found by solving a homogeneous system of linear equations. 
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Partitions 



To describe properties of v related to partitioning the system into any two subsystems, we 
need to consider all possible subsystems T and T" such that S — TUT' and the corresponding 
W and W such that V — W <SiW'. These quantities are given by 

T^Sj, T' = Si\j, W = Vj, W' = Vi\j, JeP'il), 
Sh = UhenSh, Vh = ^henVh, H e P'{I). 

Here / \ J is the relative complement of J in /. Also, P'{I) = P{I) \ {0,1}, where P{I) 
is the power set of / (the set of all subsets of /). We use P'{I) instead of P{I) to exclude 
partitions with empty subsystems {T, T') = (0, S) and (T, T') = {S, 0). 
Now, for each J e P'{I), we define the corresponding map 

fj{v) : Vj Vi\j, fj{v){w) =v®w*, 

its kernel Kj{y) = ker fj{v), and its nullity nj{v) = dim Kj{v). 

In terms of arbitrary bases {ei}i<i<dimVj and {e^}i<i<dimV/\j for the spaces Vj and Vi\j, 
we have 

dimVj dimVj\j 

Yl Yl 

1=1 j=l 

dim Vj dim V>\ j 

fjiv){w) = Y ^i^i'^i^'p 
i=l j=l 

dim Vj 

Kj{v) = e T/j: J2 ^hi'^i = 0, J e {1, . . .,dimVi\j}y 

i=l 

Example of three qubits 

For three qubits we have 

/ = {1, 2, 3}, P'{I) = {{!}, {2}, {3}, {1, 2}, {1, 3}, {2, 3}}, 
which gives the maps 

2 2 2 

/{1}(^^) -Vi^Vi^Vs, f{l}iv){w) = J21ZY1 ^h,j2,j3^he2,j2 ® 63 

jl = l 32 = 1 jz = l 

2 2 2 

f{2}{v) ■.V2^Vi® T/3, f{2}{v){w) ^YYY ^h,32,3z'^ 32^1,31 ® 63,^3, 

Jl=l j2=li3=l 
2 2 2 

31,32,33"^ 3z^'^,h ® 62,j2) 

Ji=l J2=l J3=l 
2 2 2 

/{1,2} : 14 ® 1/2 ^ 1^3, f{l,2}{v){w) = E E E ^h,32,3z'^h,32^^,3Z^ 

Jl=l J2=l J3=l 
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2 2 2 

/{1,3} ■.Vi®V^^ V2, f{l,^{v){w) = YYY^ Vj^,h,j3^juj3e2,j2, 

ii=i i2=i i3=i 
222 

/{2,3} ■.V2^V3^ Vi, f{2,3}{v){w) = ^^Yl ^h,32,3z'^32,3z^hh- 

ji=ii2=ii3=i 

Tensor products of maps 

To obtain the complete entanglement information about a vector i>, we need to describe 

its properties related to partitioning the system S' into any number of subsystems. For 
this purpose we construct the set of new maps )}jgp'(/) from the set of the maps 

{/j(^)}jgp'(/) using the operation of the tensor product. The new maps should be linear in 
V, and it should be possible to compare them with each other, for example, by comparing 
their kernels. Linearity in v requires that the only other maps allowed in the construction are 
the identity maps. Comparison of the new maps is possible only if their domains coincide, 
and a natural choice for such a common domain is the space V . These requirements fix the 
form of the new maps, 

fj{v):V ^Vi\j®Vi\j, fj{v) = fj{v)®idi\j, fj{v){x®y) = {v®x*)®y, 

where idw' '■ W — )■ W is the identity map. Jjct Kjiy) = ker fj{v) and nj{v) = dim i^'j(v) 
for each J E P'{I)- We note the relation Kj{v) — Kj{v) ® Vi\j, which follows from the 
identities 

ker {fj{v) (g) id/y) = ker fj{v) ® V7\j + Vj® ker id/y 
and kerid/\j = {0}. 

In terms of arbitrary bases {ej}i<j<dimyj and {e^}i<j<dimyj\j for the spaces Vj and Vi\j, 
we have 



i=i j=i 

dim Vj dim Vi\j dim Vi\j 

fj{v){w)^ J2 Y Yl VijWi^ke'j <^ e'k, 

i=l j=l k=l 
dim Vj 

kj{v) = eVj: ^ VijWi^k = 0, j e {1, . . . ,dimV/y}, A; e {1, . . . ,dimVf\j}|. 

i=l 

Finally, the set L{v) = {Kj{v)}j^pi(^j) is the desired set of subspaces of V that describes 
entanglement properties of v. 

Example of three qubits 

For three qubits the maps fj{v) are 

f{iy{v) : y ^2 <8) V^3 «) 1^2 (8) V^3, 
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ii=l j2=l i3=l k2 = l k-i=l 

22222 

ii=l i2 = l i3=l A:i = l ks,=l 

f{3}{v) : V ^ Vi ® V2 (g) Vi (g) V2, 
22222 

Jl=l J2=l J3=l fel=l fc2=l 

/{i,2}(v): 1^3® 1^3, 

2 2 2 2 

ii=i i2=i i3=i fc3=i 

/{i,3}(^^): V ^V2®V2, 

2222 

/(.,3,w(«') = E E E E 

'^h,32,3z'^h,k2,33^'i,32 ® ^2,k21 

jl = l j2 = l j3 = l fe2 = l 

/{2,3}M:'t"^^i<8)'t"i, 

2222 

/{2,3}(^^)(W) = 1] 1] J] ^h,32,hWk,,j2,33(^l,h ® ^hkl- 

il = l 72 = 1^3 = 1 fel = l 

Spaces and their intersections 

Results in linear algebra [22] show that the complete information about a set of linear 
subspaces is given by the dimensions of the subspaces and of all their intersections. Each 
linear space is identified by its dimension, and the intersections are needed to account for 
the relative positions of the subspaces. We specify such intersections for each set of subsets 
of 7, 

KQ{v)^nj^QKj{v), fiQiv)^ dim Kq{v), Q e P{P'{I))- 

Consequently, considering all such intersections, we find the sequence (ordered set) of new 
invariants describing all entanglement properties of v, 

N{v) = (nQ(t;))Qgp(p/(/)). 

We order elements of N{v) using the canonical ordering of elements of P{I), which is ob- 
tained from binary representations of elements of P{I) considered as {0, 1}^. We call el- 
ements of N{v) algebraic invariants of v because they are derived using standard tools of 
linear algebra. 

Finally, we define the equivalence relation ~y on the set V by setting v' ~y v if and only 
if N{v') = N{v), for each v,v' e V. The relation ~y defines the quantities C"{v), v", C", 
V" in the same manner as the relations ~y and ~y define the quantities C('i'), v, C, V and 
C'{v), v' , C", V' , respectively. 
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The proceeding development shows that the equivalence relations ~y and ~y are identical 
and proves the following theorem. 

Theorem 1. There is a one-to-one correspondence between the quotient set C and the 
sequence of values of the algebraic invariants {N{v))v^v 



Independent invariants 

In general, there are certain algebraic relations between elements of N{v). For example, 

dim Vj - - — -— = dim Vi\j - — — — 
dim Vj\j dim Vj 

is true for all v E V and J G P'{I)- We can say, for example, that nj{y) is an independent 
invariant and ni\j{v) is a dependent invariant, which can be done consistently by taking J 
only from an appropriate subset of P'{I)- 

It is convenient to remove dependent elements from N{v) by defining a subsequence of 
independent invariants 

N\v) = {hQ{v))Q^R, RCP{P'{I)), 

which we call a generating sequence of invariants of v. For consistency, we use the same 
R for each v E V. We order elements of N'[v) canonically._ All elements of N'{v) are 
algebraically independent of each other, and all elements of N{v) which are not in N'{v) 
can be algebraically expressed in terms of elements of N'(v). We choose R with the smallest 
number of elements; although this choice is not unique, all such choices are equivalent for 
our purposes. 

It remains to choose the set R. We define R = limk^oo Rk, where the sequence of sets 
{Ri, -R2, ■ ■ ■ ) is such that Rk 2 Rk+i for each k e N. We set Ri = P{P'{I)) and find the 
elements of the sequence iteratively by the following steps that remove dependent invariants: 

1. If there exist X e R'^ and Xi, G X such that Q X2, then R'^ = {R',^ \ {X}) U 
{X \ {^1}) for any such X, Xi, otherwise, R'^ — R'^. 

2. If there exists X e R'l such that XiHXs = for any Xi, X2 G X, then i?'^" = Rk\{X} 
for any such X; otherwise, R'^' — R'^. 

3. If there exist X,Y e R';^' such that Xi e Y for any Xi e X, then Rk+i = R^' \ {X} 
for any such X; otherwise, Rk+i — R'^'- 

If there is more than one choice for X (and for Xi in step 1) that satisfies the conditions in a 
given step, then any such choice can be made. (The resulting sequence {Ri, R2, ■ ■ ■) depends 
on these choices.) For any such choice, however, the sequence is convergent and its limit 
R = limfc^oo Rk is reached after a finite number of iterations, i.e. there exists m G N such 
that Rk = R for each k > m. Even though the above choices can lead to different sets R 
and resulting generating sequences N'{v), they lead to the same entanglement classification. 
This completes the construction of each generating sequence of invariants N'{v). 
The relation 

Kq{v) = njeQ{Kj{v) Vr\j), Q G P{P\I)) 
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implies 



nQ{v)^nQ{v)dimVj\uj^Qj, Uj^qJ^I, Q e P{P'{I)), 

where we set dim V0 = 1. (To prove this, note that Kq(v) does not involve Kj'[v), where J' G 
/\UjeQj.) For n >3, this relation between the invariants means that such ngiv) describes 
properties of v related to partitioning the system into at most |/\ Uj^q J\ subsystems. For 
such cases, it is convenient to replace 'hQ{v) with nqiv) and define the set of invariants 



nqiv), Uj^qJ C I, 



N"{v) = {mQ{v))Q^R, mgiv) 
We order elements of N"{v) canonically. We give our explicit solutions in terms of N"{v). 
General forms of states 

As our main computational device, we use the general forms of elements of V". We obtain 
them from expressions for elements of ker /(v) for a map f(v) : W W, to derivation of 
which we now turn. We choose arbitrary bases {'Ui}i<j<dimM/ and {'Ui}i<i<dim w for the 
spaces W and W, respectively, and represent a vector v G V in terms of its coordinates, 

dim TV dim W 

J2 Yl ^^'i^^ ® ^'p "L^^^j} C F. 
i=i j=i 

It follows that V decomposes according to 

dimW dimW 

i=l j=l 

dimW dimW 
V — ^ Uj ® u'j, Uj — ^ VijUi, {Uj} C W. 
j=l i=l 

The defining relation v <S> w* — ior w E kei f{v), which is a system of homogeneous linear 
equations for the coordinates of w, now implies the general form of v, 

dim W—n{v) 

V = ^ (g) {wi} C W, {w-} C W\ 

i=l 

dim span {{wi}) — dim span {{w'^}) — dimH^ — n{v), 

where n{v) — dimker /(f ) and the dimension of the span of a set of vectors is the number of 
its linearly independent elements. This decomposition is unique up to linear transformations 
Wi I— 7- J2j B'ijWj and w[ 1— )> J^j^ij'^p where B and B' are nonsingular square matrices of 
order dim 14^ — n{v) that satisfy the condition B^B' = 1. 

When considering the above general forms of elements of V resulting from different choices 
of W and W such that V — W <Si W, we need to choose {wi} and {w-} (using appropriate 
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B and B') such that the corresponding decompositions are consistent for all such choices. 
This results in restrictions on allowed values of the invariants in N{v) and, consequently, 
leads to the classification of all entangled states. 

The described method solves the restricted classification problem for arbitrary 
Obtaining explicit solutions, however, is entirely different matter. We did not obtain such 
solutions for arbitrary {Vi}jg/, but we found them for numerous examples given in the 
following section. 

Particularly interesting are cases where the spaces in {Viji^i are of equal dimensions. 
The resulting permutation symmetry among the spaces reduces the equivalence classes to 
sets of classes related by the symmetry. As a result, representative elements for the sets of 
classes take simple forms. We have explicit solutions for two such symmetric examples. 

IV. EXAMPLES 

The present classification method works for arbitrary finite n and D — (dimVi)jg/. The 
case n = 2 is easily solved [21] for arbitrary D. We now apply our method to the case n — 3 
for a large selection of values of D and the case n — A, D — {2, 2, 2, 2) (four qubits). 

A. n — 3 

Independent invariants for n = 3 are given by the sets 

Qi = {{l}}, Q2 = {{2}}, g3 = {{3}}, g4 = {{l,2},{l,3},{2,3}}. 

The sets Qi,Q2,Q3 and Q4 lead to invariants related to partitioning the system into two 
and three subsystems, respectively. For each of these invariants, there are corresponding 
invariants generated by the transpose maps, which do not need to be considered. Since all 
other partitions lead to dependent invariants, we choose the generating set of invariants 

N"{v) = {nQ,{v),nQ^{v),nQ^{v),nQ^{v)) 

for each v eV. 

For the set of equivalent classes we find 

C" = {Co}U{Ck,,k,,k3j- ^1 e {l,---,di},k2 e {I,...,d2},k3 G {l,...,d3},j e Mk,,k2,ks}^ 

where D = {di, d2, ds) and Mk^^k2M a certain set of natural numbers that is symmetric in 
^1) ^3- The values of the invariants in N"{v) for the classes Cq and Ckj^.k-zM-j given in 
Table I. Although we do not have a general formula for ^j.fea fo^ arbitrary (fci, k2, k^), we 



TABLE I. The values of the invariants in N"{v) for n = 3, D = (^1,^2, d^). 













Co 


di 


d2 


ds 


did2d3 


Cki,k2,k3,j 


di — ki 


d2 - k2 


ds - ks 


did2d3 - kidi - k2d2 - ksds + {Mk^^k^M)] 



give Mk^^k2M various particular values of (/ci, A;2, in Table II, which is our main result 
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for the case n = 3. With analogous computations for additional values of {k1.k2.k3), the 
table can be easily expanded. Such a table is directly used for explicit computations of C" 
for various values of D. In particular, the values of Mk^^k2,kz given in Table II suffice to find 
the set of classes C" for each value of D given in Table III; the latter table gives only the 
number of classes \C"\. As illustrative examples and because of space limits, wc present here 
the full results only for D = {2,2, d) and D = {2,3, d), where d is arbitrary, in Tables IV 
and V, respectively. For the symmetric case D = (2, 2, 2), there are 5 sets of classes related 
by permutations of {Vi, V2, V3}; Table VI lists the sets and their representative elements. 



TABLE II. The set Mfc^^^jjiks for various values of {ki, k2, ks). The notation (m, . . . , m') means all 
integers between and including m and m'. 







{kx,k2,kz) 




1,1,1) 


(2) 


(2,6,6) 


(7,..., 23, 28, 29) 


1,2,2) 


(5) 


(2,6,7) 


(5,..., 22, 24, 25, 26, 34) 


1,3,3) 


(10) 


(2,6,8) 


(8,..., 25, 31) 


1,4,4) 


(17) 


(2,6,9) 


(13, 14, 16, ... , 20, 22, ... , 26, 29, 30) 


1,5,5) 


(26) 


(2,6, 10) 


(20,23,24,28,29,31) 


!.(),()) 


(37) 


(2,6,11) 


(29,34) 


2,2,2) 


(4,5) 


(2,6,12) 


(40) 


2,2,3) 


(5,6) 


(3,3,3) 


(2,. ..,8,10) 


2,2,4) 


(8) 


(3,3,4) 


(2,..., 11) 


2,3,3) 


(4,. ..,8) 


(3,3,5) 


(2,..., 11, 14) 


2,3,4) 


(5,..., 8, 10) 


(3,3,6) 


(2,...,12) 


2,3,5) 


(8,10) 


(3, 3, 7) 


(4,..., 11, 14) 


2,3,6) 


(13) 


(3,3,8) 


(10,11,14) 


2,4,4) 


(5,..., 13) 


(3,3,9) 


(18) 


2,4,5) 


(5,..., 13, 16) 


(3,4,4) 


(2,..., 12, 14) 


2,4,6) 


(8,9,10, 12,13, 15) 


(3,4,5) 


(2,. ..,16) 


2,4,7) 


(13,16) 


(3,4,6) 


(2,..., 16, 20) 


2,4,8) 


(20) 


(3,4,7) 


(2.. ..,17) 


2,5,5) 


(6,..., 16, 19, 20) 


(3,4,8) 


(2,..., 17, 19) 


2,5,6) 


(5,..., 18, 24) 


(3,4,9) 


(2,..., 17, 20) 


2,5,7) 


(8, 10,..., 18, 20, 22) 


(3,4, 10) 


(5,..., 16, 19,20) 


2,5,8) 


(13,15,16, 19,20,22) 


(3,4,11) 


(14, 16,20) 


2,5,9) 


(20, 24) 


(3,4, 12) 


(25) 


2,5,10) 


(29) 







It is easy to obtain general expressions for M^^^^^M for various particular values of 
{ki, k2, k^), and we give here just a few such results: 

Mk,MMk2-i ^{...,kl + kl- 2{ki + k2) + 5,kl + kl- {ki + k2) + 2). 

These and similar readily available expressions for M/jj ^2,^3 suggest certain patterns, which 
might eventually lead to the general result for arbitrary {ki, k2, k^). 
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TABLE III. The numbers of equivalence classes 


\C"\ for n 


= 3 and various values of D. 




D 




\C"\ 


D 


ic 1 


U 


|C 


(2,2,2) 




7 


(2,5,5) 


7^ 

11 


/o o o ^ 
(o, o, oj 


on 
39 


(2,2,3) 




9 


(2,5,6) 


yy 




DU 


(2,2,d) 


d>4 


10 


(2,5,7) 


1 1 o 


/o o c ^ 
(o, oj 


75 


(2,3,3) 




17 


(2,5,8) 


1 OA 




00 
00 


(2,3,4) 




23 


(2,5,9) 


izz 


(3,3, 7) 


97 


(2,3,5) 




25 


(2,5,d),d> 10 


1 OQ 




iUU 


(2.3,r/) 


d > G 


2G 


(2,6,6) 


1/11 


(o. o, , d \) 


11) 1 


(2,4,4) 




39 


(2,6,7) 


1 TT 

177 


(3,4,4) 


1 no 
103 


(2,4,5) 




51 


(2,6,8) 


onQ 
/Uo 


4, Oj 




(2,4,6) 




58 


(2,6,9) 


oi n 
219 


(3,4,6) 


1 'TO 

17o 


(2,4, 7) 




60 


(2,6,10) 




(3,4, 7) 


one: 


(2,4,rf) 


,d>8 


61 


(2,6,11) 


228 


(3,4,8) 


226 








{2,6, d),d> 12 


229 


(3,4,9) 


244 












(3,4,10) 


258 












(3,4,11) 


261 












(3,4,d),d> 12 


262 



TABLE IV. The entanglement classes, their algebraic invariants, and their representative elements 
for = 3, D = {2,2, d). Classes for which any of the invariants in the set N"{v) are negative 
should be discarded. Classes within a horizontal block are added each time d increases by 1, so 
that there are 7, 9, 10 classes for d = 2, d = 3, d > 4, respectively. Each expression [71,^2,^3] 
stands for uij^ (8) U2,j2 ® ^43,^3, where {uij} is a set of any linearly independent elements of V^. 



X"{r) 



Co 


{2,2,d,Ad) 





Ci 


(l,l,d-l,3a!-2) 


[1,1,1] 


C2 


(0,0,d- l,3(i- 3) 


[1,1,1] + [2,2,1] 


C3 


(0, l,d-2,2ci-l) 


[1,1,1] + [2, 1,2] 


Ca 


(l,0,(i-2,2d- 1) 


[1,1,1] + [1,2, 2] 


Cn 


(0,0,d-2,2d-3) 


[1,1,1] + [1,2, 2] + [2, 1,2] 


Cf, 


(0,0,d-2,2d-4) 


[1,1,1] + [2, 2, 2] 


C7 


(0,0,d-3,(i-2) 


[1,1,1] + [1,2,2] + [2,2,3] 




(0,0,d-3,d-3) 


[1,1,1] + [1,2,2] + [2,1,2] + [2,2,3] 


C9 


(0,0,d-4,0) 


[1,1,1] + [1,2, 2] + [2, 1,3] + [2, 2, 4] 



The needed computations for the above cases are lengthy but elementary, and we do 
not give their details here. Instead, we invite the reader to study graphical representation 
of entanglement classes for the cases D — (2, 2, d) and D — (2, 3, d) in Figs. 1 and 2, 
respectively, which can be easily generalized for arbitrary D. Although these and similar 
figures cannot replace the actual computations, they are useful in understanding relations 
between the classes, finding their general properties, and, perhaps, even in solving the general 
case. In this regard, generalizations of Table VI seems to be particularly promising when 
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TABLE V. The entanglement classes, their algebraic invariants, and their representative elements 
for n = 3, = (2, 3, d). Classes for which any of the invariants in the set N"{v) are negative 

should be discarded. Classes within a horizontal block arc added each time d increases by 1, so 
that there are 9, 17, 23, 25, 26 classes for d = 2, d = 3, d = 4, d = 5, d > 6, respectively. Each 
expression [ji,i2, js] stands for uij^ ®'"2j2 ®^3j3) where {uij} is a set of any linearly independent 
elements of Vi. 











V 














1 <Jj U/j Ul* J 



















(1-, 


d 2/7—1 


■id 






]_] 














Oil 






1 " 
1 


1 




1 1 
1 








ArJ 






1 " 
1 


1 


\0 1 


91 




r\ 
04 




Arl 






1 
1 


1 


fl 9 


91 






1 // — 9 


Arl 

"-til 






1 " 
L 




fl 9 


91 4- f9 1 91 




(la 




Ad 






]_] 


_|_ 


[2 2 


2] 






rn r7 — 9 


Ad 






1 " 
L 


_l_ 


fl 9 


91 -1- f9 11 
1 ^ 1 , '-^ , 1 






i{\ W — 9 


Ad 






1 " 
± 


_l_ 
\ 


fl 9 

, 


91 -1- f9 9 11 


4- f9 91 




now 


OIL 


1 ^ 


fl 1 


1 " 
1 


1 


fl 9 
, 


91 4- fl "^l 

z| 11,0,01 








Oil 


— A\ 


fl 1 


1 1 
± 


1 


fl 9 

f , 


91 4- f9 1 "^l 








OUi 




fl 1 


11 
± 


_|_ 


fl 9 


91 4- f9 1 91 


-1- f9 9 "^1 




(0, 0, d — 3, 


3d 




fl 1 




-1- 


fl 2 


91 4- fl Si 


+ f2 1 21 




(0, 0, d — 3 


3d 




fl 1 


1" 


\ 


fl 2 


91 4- f9 "^l 




Ci4 


(0,0,d- 3 


3d 


-7) 


[1,1 


1; 


+ 


[1,2 


2] + [1,3, 3] 


+ [2, 1,2] + [2, 2, 3] 


Cl5 


(0,0,(i- 3, 


3d 


-8) 


[1,1 


1] 


+ 


[1,2 


2] + [2,1,3] 


+ [2,3,1] 


C16 


(0,0,d-3, 


3d 


-9) 


[1,1 


1; 


+ 


[1,2, 


2] + [2, 2, 2] 


+ [2,3,3] 


Ci7 


(0,l,d-4, 


2d 


-2) 


[1,1 


1] 


+ 


[1,2 


2] + [2, 1,3] 


+ [2,2,4] 


C18 


(0,0,d-4 


2d 


-3) 


[1,1 


1; 


+ 


[1,2 


2] + [1,3,3] 


+ [2,3,4] 


C19 


(0,0,(i-4 


2d 


-5) 


[1,1 


1; 


+ 


[1,2 


2] + [1,3,3] 


+ [2,2,4] + [2,3,1] 


^^20 


(0,0,(i-4 


2d 


-6) 


[1,1 


1; 


+ 


[1,2 


2] + [2,2,3] 


+ [2,3,4] 


C2I 


(0,0,d-4, 


2d 


-7) 


[1,1 


1; 


+ 


[1,2, 


2] + [1,3, 3] 


+ [2, 2, 3] + [2, 3, 4] 


C22 


(0,0,d-4, 


2d 


-8) 


[1,1 


1] 


+ 


[1,2 


2] + [1,3, 3] 


+ [2, 1,2] + [2, 2, 3] + [2, 3, 4] 


(^23 


(0,0,d- 5, 


d- 


3) 


[1,1 


1; 


+ 


[1,2 


2] + [1,3,3] 


+ [2, 1,4] + [2, 2, 5] 


C24 


(0,0,d- 5, 


d- 


5) 


[1,1 


1; 


+ 


[1,2 


2] + [1,3,3] 


+ [2, 1,3] + [2, 2, 4] + [2, 3, 5] 


(^25 


(0,0,d-6, 


0) 




[1,1 


1; 


+ 


[1,2 


2] + [1,3,3] 


+ [2, 1,4] + [2, 2, 5] + [2, 3, 6] 



solving the symmetric case D ~ {d, d, d) for arbitrary d. 
B. n = 4 

Table VII lists sets that give independent invariants for n — A, arranged according to 
types of partitions of the system. The sets Qi, . . . , (^4 and Q5, . . . , Qs lead to invariants 

related to partitioning the system into two and three subsystems, respectively. Partitions 
into four subsystems are of three different types and are given by the sets Qg, . . . , (5i4, 
and Q15, • • -Qis, and Qig- For each of these five types, there are corresponding invariants 
generated by the transpose maps, which do not need to be considered. Since all other 
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TABLE VI. Representative elements for the sets of equivalence classes for n = 3, D = (2, 2, 2) 
induced by the permutation symmetry of the spaces in {Vi,V2)^}- A representative element v 
is given hy v = Avi, where A: V ^ V is a certain linear operator and f i G y is a fixed vector. 
(Without loss of generality and for comparison with other tables, we choose vi = [1,1,1].) The 
operator Oj is defined by ai[. . . , 1, . . . ] = [. . . , 2, . . . ] and ai[. . . ,2, . . .] = [. . . , 1, . . . ], where only 
the ith index changes. To obtain all classes in each group, all possible choices of the indices 
{i,j,k} = {1,2,3} should be considered. 

A 

'Co 
Ci 1 

{C2,C3,C4} l + aiGj 

C5 1 + ai{aj + Gk) 

Ce 1 + a^ajak 



i 

C3 



C9 



FIG. 1. Graphical representation of entanglement classes for n = 3, D = (2,2, d). Each vertex 
corresponds to a certain expression [ji,j2,j3] in a representative element v for each class, and v 
is the sum of such expressions over all vertices of a given three-dimensional lattice; see Table IV. 
The invariants nQ^(u), nQ^{v), nQ^{v) equal the numbers of linearly independent two-dimensional 
lattices. To obtain the corresponding representations for D = (2,2,d), we remove 4 — d cubes 
from tops of stacks for d < 4 or add d — 4 cubes on top of stacks without adding any new vertices 
for d > 4. There are only 10 classes for any d > 4 because the number of linearly independent 
two-dimensional lattices along one of the directions is already maximal (four) for the class Cg. The 
construction for arbitrary D is analogous; see, for example. Fig. 2. 



partitions lead to dependent invariants, we choose the generating set of invariants 
N"{v) = {nQ^{v), . . . ,nQ^{v),hQ^{v), . . . ,hQ,,{v)) 

for each v & V. 

As an illustrative example, we take D = (2,2,2,2). It is convenient to partition C" into 
three sets, 

c" = c'l, u c'l, u c;,3, 

according to possible forms of representing elements for classes in each set. The set Cji 
consists of classes that can be represented by elements with coefficients in linear combinations 
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Co 
to 
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to 

Ci 
to 

Cio 
to 

C19 



to 

to 

to 
C20 



to 

C3 
to 

C12 
to 



to 

C4 
to 

to 



to 

C5 
to 

to 



to 

to 

to 

11^ 



to 

11^ 

to 

C16 

to 

II- 



to 

m 

Cs 
to 

11^ 



C 



17 



G 



21 



22 



G 



23 



c. 



24 



25 



FIG. 2. Graphical representation of entanglement classes for n 
further details. 



3, D = (2,3,d). See Fig. 1 for 



TABLE VII. Sets that give independent invariants for n = 4. 



Qi = {{1}} Q9 = {{1,2},{1,3,4},{2,3,4}} 

Q2 = {{2}} Qio = {{1,3},{1,2,4},{2,3,4}} 

0,3 = {{3}} On = {{1,4},{1,2,3},{2,3,4}} 

O4 = {{4}} O12 = {{2,3},{1,2,4},{1,3,4}} 



O5 = {{1,2}, {1,3}, {2, 3}} Oi3 = {{2,4},{1,2,3},{1,3,4}} 

Qe = {{1,2},{1,4},{2,4}} Q14 - {{3, 4}, {1, 2, 3}, {1, 2, 4}} 



O7 = {{1,3},{1,4},{3,4}} Qi5 = {{1,2},{1,3},{1,4},{2,3,4}} 

Os = {{2,3},{2,4},{3,4}} Q16 = {{1,2},{2,3},{2,4},{1,3,4}} 

Oi7-{{l,3},{2,3},{3,4},{l,2,4}} 
Oi8 = {{l,4},{2,4},{3,4},{l,2,3}} 



0i9 - {{1,2,3},{1,2,4},{1,3,4},{2,3,4}} 



of bases vectors taken from {0, 1}. The set Cj 2 consists of classes that do not belong to Cj\ 
and that can be represented by elements with coefficients in linear combinations of bases 
vectors taken from {0, 1, —1}. The set Cj^ consists of classes that do not belong to either 
C'j' I or Cj 2- The classes in Cji are the simplest and the most typical, and the classes in Cj^ 
are the most complex and the least typical. It is clear that classes in Cj\ and Cj2 can be 
represented by elements with coefficients in linear combinations of bases vectors taken from 
other sets besides {0, 1} and other sets besides {0, 1}, {0, 1, —1}, respectively. Nevertheless, 
our results show that the chosen partition of C" is by no means arbitrary. 
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We find 



~ {^0; • • • ; C29, C34, . . . , Cee, Ces, • • • , C'82}, Cj 2 — {C'so; C'si; C'32, e'er}; C'/'s 5 {C'ss}- 

We used the Monte Carlo method to search for the set C" ^, and it is possible that it contains 
additional classes besides C33. However, our results show that such additional classes are 
very rare with respect to a measure that is uniform on the space of coefficients in linear 
combinations of bases vectors. Tables VIII, IX, X fist the classes, their independent invari- 
ants and representative elements. In these tables, all 83 classes appear in 27 fundamental 
sets of classes related by permutations of {Vi, V2, V3, V4}. Table XI lists the sets of classes 
and their representative elements. 

TABLE VIII. The entanglement classes, their independent algebraic invariants, and their rep- 
resentative elements for n = 4, D = (2,2,2,2). Each expression [ii, J2,i3, J4] stands for 
uiji <^'W2,j2 '^'"3j3 ® ^4,j4) where {uij} is a set of any linearly independent elements of Vi. Classes 
within each horizontal block are related by a permutation symmetry of {V^jig/. 



N"{v) V 



Co 


(2 


2, 


2 


2, 


8 


8 
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8, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16) 
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1, 


1 


1, 


4 


4 


4 


4,10,10,10,10,10,10,8,8,8,8,11) 


[1, 


1 


1, 
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C2 


(0 


0, 


1 


1, 


3 


3 


2 


2,9,7, 7,7,7,10,3,3,7,7,10) 


[1, 


1 


1, 


1 


+ 


[2 


2 


1 


1; 




C3 


(0 


1, 





1, 


3 


2 


3 


2,7,9,7,7,10,7,3,7,3,7,10) 


[1, 


1 


1, 


1; 


+ 


[2 


1 


2 


1; 




Ci 


(0, 


1, 


1, 


0, 


2, 


3 


3, 


2,7,7,9, 10,7,7,3,7,7,3, 10) 


[1, 


1, 


1, 


1; 


+ 


[2 


1, 


1 


2; 




C5 


(1 


0, 





1, 


3 


2 


2 


3,7,7, 10,9,7,7,7,3,3,7,10) 


[1, 


1 


1, 


1 


+ 


[1 


2 


2 


1; 




Ce 


(1 


0, 


1 


0, 


2 


3 


2 


3,7,10,7, 7,9,7, 7,3,7,3,10) 


[1, 


1 


1, 


1 


+ 


[1 


2 


1 


2; 




Cv 


(1 


1, 





0, 


2 


2 


3 


3,10,7, 7, 7,7,9,7, 7,3,3,10) 


[1, 


1 


1, 


1; 


+ 


[1 


1 


2 


2; 




Cs 


(0 


0, 





0, 











0,0,0,9,9,0,0,0,0,0,0,9) 


[1, 


1 


1, 


1; 


+ 


[1 


2 


2 


1; 


+ [2, 1,1, 2] + [2, 2, 2, 2] 


C9 


(0 


0, 


0, 


0, 


0, 





0, 


0,0,9,0,0,9,0,0,0,0,0,9) 


[1, 


1, 


1, 


1; 


+ 


[1 


2, 


1 


2; 


+ [2, 1,2,1] + [2, 2, 2, 2] 


Cia 


(0 


0, 





0, 











0,9.0,0,0,0,9,0,0,0,0,9) 


[1, 


1 


1, 


1; 


+ 


[1 


1 


2 


2; 


+ [2, 2, 1.1] + [2, 2, 2, 2] 


Cn 


(0 


0, 





1, 


1 


2 


2 


2,5,5,7,5,7, 7,2,2,2,7, 8) 


[1, 


1 


1, 


1 


+ 


[2 


1 


2 


1 


+ [2,2,1,1] 


C12 


(0 


0, 


1 


0, 


2 


1 


2 


2,5,7,5,7,5,7,2,2,7,2,8) 


[1, 


1 


1, 


1; 


+ 


[1 


2 


1 


2; 


+ [2,2,1,1] 


Ci3 


(0 


1, 





0, 


2 


2 


1 


2,7,5,5,7,7,5,2,7,2,2,8) 


[1, 


1 


1, 


1; 


+ 


[1 


1 


2 


2; 


+ [2,1,2,1] 


(7i4 


(1, 


0, 


0, 


0, 


2, 


2 


2, 


1,7,7,7,5,5,5,7,2,2,2,8) 


[1, 


1, 


1, 


1; 


+ 


[1 


1, 


2 


2; 


+ [1,2,1,2] 


Ci5 


(0 


0, 


0, 


1, 


0, 


2 


2, 


2,4,4,7,4,7,7,2,2,2,7,7) 


[1, 


1 


1, 


1; 


+ 


[2 


2 


2 


1; 




C16 


(0 


0, 


1 


0, 


2 





2 


2,4,7,4,7,4,7,2,2,7,2,7) 


[1, 


1 


1, 


1; 


+ 


[2 


2 


1 


2; 




Ci7 


(0 


1, 





0, 


2 


2 





2,7,4,4,7,7,4,2,7,2,2,7) 


[1, 


1 


1, 


1 


+ 


[2 


1 


2 


2; 




Cis 


(1 


0, 





0, 


2 


2 


2 


0,7,7, 7,4,4,4, 7,2,2,2,7) 


[1, 


1 


1, 


1; 


+ 


[1 


2 


2 


2; 




Cl9 


(0 


0, 


0, 


0, 


1, 


1 


1, 


1,5,5,5,5,5,5,2,2,2,2, 7) 


[1, 


1, 


1, 


1; 


+ 


[2 


1, 


1 


2; 


+ [2, 1,2,1] + [2, 2, 1,1] 


C20 


(0 


0, 


0, 


0, 


1, 





0, 


1,2,2,4,5,2,2,1,1,1,1,6) 


[1, 


1, 


1, 


1; 


+ 


[1 


2, 


2 


1; 


+ [2,2,1,2] 


C21 


(0 


0, 





0, 





1 


1 


0,2.2,5,4,2,2,1,1,1,1,6) 


[1, 


1 


1, 


1 


+ 


[2 


1 


1 


2; 


+ [2,2,2,1] 


C22 


(0 


0, 





0, 





1 





1,2,4,2,2,5,2,1,1,1,1,6) 


[1, 


1 


1, 


1; 


+ 


[1 


2 


1 


2; 


+ [2,2,2,1] 


C23 


(0 


0, 





0, 


1 





1 


0,2,5,2,2,4,2,1,1,1,1,6) 


[1, 


1 


1, 


1; 


+ 


[2 


1 


2 


1; 


+ [2,2,1,2] 


C24 


(0, 


0, 


0, 


0, 


0, 





1, 


1,4,2,2,2,2,5,1,1,1,1,6) 


[1, 


1, 


1, 


1; 


+ 


[1 


1, 


2 


2; 


+ [2,2,2,1] 


C25 


(0, 


0, 


0, 


0, 


1, 


1 


0, 


0,5,2,2,2,2,4,1,1,1,1,6) 


[1, 


1, 


1, 


1; 


+ 


[2 


1, 


2 


2; 


+ [2,2,1,1] 


C26 


(0 


0, 


0, 


0, 


0, 





0, 


0,4,4,4,4,4,4,0,0,0,0,6) 


[1, 


1, 


1, 


1; 


+ 


[2 


2, 


2 


2; 
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TABLE IX. The entanglement classes, their independent algebraic invariants, and their rep- 
resentative elements for n = 4, D = (2,2,2,2). Each expression [^1,^2,^3,^4] stands for 
<^ ■"2,^2 ® ^3,i3 ""4^4 5 where {uij} is a set of any linearly independent elements of Vi. Classes 
within each horizontal block are related by a permutation symmetry of {V^jjg/. 



N"{v) V 


C27 (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 5, 0, 0, 0, 0, 0, 0, 6) [1, 1, 1, 1] + [1, 1, 2, 2] - 
C28 (0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 5, 0, 0, 0, 0, 0, 6) [1, 1, 1, 1] + [1, 1, 2, 2] - 
C29 (0, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 5, 0, 0, 0, 0, 6) [1, 1, 1, 1] + [1, 2, 1, 2] - 


f [1,2, 1,2] H 
f [1,2, 2,1] H 
f [1,2, 2,1] H 


h[2,l,2,l]H 
h[2,l,l,2]H 
h[2,l,l,2]H 


h [2, 2, 1,1] 
h [2, 2, 1,1] 
h [2, 1,2,1] 



C30 (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 6) 



C31 (0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 6) 



C32 (0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 6) 



[1, 1, 1, 1] - [1, 1, 1, 2] - [1, 1, 2, 1] - [1, 2, 1, 1] + [1, 2, 2, 1] 
+ [1, 2, 2, 2] - [2, 1, 1, 1] + [2, 1, 1, 2] + [2, 1, 2, 2] + [2, 2, 1, 2] 

+ [2, 2, 2,1] + [2, 2, 2, 2] 

[1, 1, 1, 1] - [1, 1, 1, 2] - [1, 1, 2, 1] - [1, 2, 1, 1] + [1, 2, 1, 2] 
+ [1, 2, 2, 2] - [2, 1, 1, 1] + [2, 1, 2, 1] + [2, 1, 2, 2] + [2, 2, 1, 2] 
+ [2, 2, 2,1] + [2, 2, 2, 2] 

[1, 1, 1, 1] - [1, 1, 1, 2] - [1, 1, 2, 1] + [1, 1, 2, 2] - [1, 2, 1, 1] 
+ [1, 2, 2, 2] - [2, 1, 1, 1] + [2, 1, 2, 2] + [2, 2, 1, 1] + [2, 2, 1, 2] 
+ [2, 2, 2,1] + [2, 2, 2, 2] 



C33 (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6) 



[1, 1,1,1] + c[l, 1,2,2] - 
+c[2,2,l,l] + [2,2,2,2]. 



(1 + c)[l, 2,1,2] - 
cei^, c^{-2,■ 



-(l+c)[2,l,2,l] 
-1,0,1} 



C34 (0, 

C35 (0, 

C36 (0, 

C37 (0, 



0,0,0,0 

0,0,0,0 
0,0,0,0 
0,0,0,1 



0,0 

0,1 
1,0 
0,0 



1,2 

0,2 
0,2 
0,2 



2,2,2,2,2,0 

2,2,2,2,2,1 
2,2,2,2,2,1 

^7 ^7 ^ 



1,1,1,5) 

0,1,1,5) 
1,0,1,5) 
1,1,0,5) 



[1,1,1,1] + 

[1,1,1,1] + 
[1,1,1,1] + 
[1,1,1,1] + 



1,2 

1,1 
1,1 
1,1 



+ 



2,1,1,1] + 

1,2,1,1] + 
1,2,2,2] + 
2,1,2,2] + 



[2,2,1, 

[2,2,2, 
[2,1,1, 
[2,2,1, 



^38 (0, 
C39 (0, 

C40 (0, 



0,0,0,0 
0,0,0,0 

0,0,0,0 



0,0 
0,0 

0,0 



0,2 
0,2 

0,4 



2,4,4,2,2,0 
4,2,2,4,2,0 

2,2,2,2,4,0 



0,0,0,5) 
0,0,0,5) 

0,0,0,5) 



[1,1,1,1] + 
[1,1,1,1] + 

[1,1,1,1] + 



1,2 
1,2 
1,1 



+ 



2,1,1,2] 
2,1,2,1] 

2,2,1,1] 



C41 (0, 

Ci2 (0, 
C43 
C44 

C45 (0, 

C46 (0, 



(0: 

(0. 



0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 

0,0,0,0 
0,0,0,0 



0,0 
0,0 
0,0 
0,0 
0,0 
0,0 



0,0 
0,0 
0,0 
0,0 

0,4 
0,5 



0,4,5,0,0,0 
0,5,4,0,0,0 
4,0,0,5,0,0 
5,0,0,4,0,0 

0,0,0,0,5,0 
0,0,0,0,4,0 



0,0,0,5) 
0,0,0,5) 
0,0,0,5) 
0,0,0,5) 
0,0,0,5) 
0,0,0,5) 



[1,1,1,1] + 
[1,1,1,1] + 
[1,1,1,1] + 
[1,1,1,1] + 

[1,1,1,1] + 
[1,1,1,1] + 



1, 1 
1,1 
1,1 
1,1 
1,1 
1,2 



+ 



1,2,1,1] + 
1,2,2,2] + 
1,2,1,1] + 
1,2,2,2] + 
1,1,2,1] + 
1,2,2,2] + 



2,1,2, 
2,1,1, 
2,1,2, 
2,1,1, 

2,2,1, 
2,1,1, 



2] + [2, 2, 1,2] 
1] + [2, 2, 2,1] 
2] + [2, 2, 2,1] 
1] + [2, 2, 1,2] 
2] + [2, 2, 2,1] 
1] + [2, 1,2, 2] 



Ci7 (0, 
C48 (0, 
Ci9 (0, 



0,0,0,0 
0,0,0,0 
().(). 0,0 



0,0 
0,0 

u.o 



0,0 
0,0 

0,4 



0,4,4,0,0,0 
4,0,0,4,0,0 
0.0.0,0,4.0 



0,0,0,5) 
0,0,0,5) 

0.0, ().:)) 



[1,1,1,1] + 
[1.1,1.1] + 
[1.1,1.1] + 



1,2 

1-2 
1.2 



2,1,1,2]- 

2,1.2,1] 

2,1.2.2: 



2,2,2, 
2,2.1, 
2,2. 1, 



C50 (0, 

C51 (0, 

C52 (0, 

C53 (0, 



0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 



0,0 
0,0 
0,0 
0,0 



0,1 
0,1 
0,2 
0,2 



1,2,1,2,2,0 
2,1,2,1,2,0 
1,1,2,2,1,0 
2,2,1,1,1,1 



0,0,1,4) 

0,1,0,4) 
1,0,0,4) 
0,0,0,4) 



[1,1,1,1] + 

[1,1,1,1] + 
[1,1,1,1] + 
[1,1,1,1] + 



1,2 
1,2 
1,1 
1,1 



2,1,2,1] 
2,1,2,2]- 
2,1,2,1]- 
1,2,1,2]- 



2,2,1, 
2,2,1, 
2,2,1, 
2, 2, 2, 



C54 (0, 

C55 (0, 

C56 (0, 

C57 (0, 

C5S (0, 

C59 (0, 



0,0,0,0 
0,0,0,0 

0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 



0,0 
0,0 

0,0 
0,0 
0,0 
0,0 



0,1 
0,2 

0,2 
0,2 
0,2 
0,2 



2,2,2,2,2,0 
1,2,2,2,2,0 

2,1.2,2,2,0 
2,2, 1,2,2,0 
2,2,2,1,2,0 
2,2,2,2,1,0 



0,0,0,4) 
0,0,0,4) 

0,0,0,4) 
0,0,0,4) 
0,0,0,4) 
0,0,0,4) 



[1,1,1,1] + 
[1,1,1,1] + 

[1,1,1,1] + 
[1,1,1,1] + 
[1,1,1,1] + 
[1,1,1,1] + 



1,2 

1,1 
1,1 
1,1 
1,1 
1,1 



+ 



+ 



2,1,2,1] + 
2,1,2,1] + 

2,1,1,2] + 
1,2,2,1] + 
1,2,1,2] + 
1,2,1,2] + 



2,2,1,1] 

1] 
1] 
1] 
1] 
1] 



2,2,1, 
2,2.1, 
2,2,1, 
2,2,1, 
2,1,2, 
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TABLE X. The entanglement classes, their independent algebraic invariants, and their rep- 
resentative elements for n = 4, D = (2,2,2,2). Each expression [^1,^2,^3,^4] stands for 
® ■"Sj^ ® ^3J3 ""4^4 5 where {uij} is a set of any linearly independent elements of Vi. Classes 
within each horizontal block are related by a permutation symmetry of {V^jig/. 







V 




















[111 
1 ± , 1 , ± 


1 
1 


^[19 
+ [i, Z 


1 9 


4- [1 


9 9 1 
Z, Z, 1 


4- f9 1 1 9 

+ [z, r, -L, z 


4- f9 1 9 11 

+ [z, r, z, ij 






_i_ro 

+ ['^, ^■ 




z, 


91 














(^U,U,U,U,U,U,U,U,U,U,Z,Z,U,U,U, U, U,U,4j 




1 
i 


+ [1, 1 


9 9 

z, z 


4_ n 
+ [1 


9 1 1 

z, J-, J-, 


4_ [9 1 1 1 

+ [z, 1, i, 1, 


4_ [9 1 9 11 

+ [z, 1, z, iJ 






_L To 9 


1 

i, 


91 












(-^62 


^u, u, u, u, u, u, u, u, u, z, u, u, z, u, u, u, U, U, 4J 


ri 1 1 


1 
1 


_L n 1 
+ [1, 1 


1 9' 


4_ fl 
+ [1 


9 9 1' 

z,z, 1 


4_ r9 1 1 1 ' 

+ [z, 1, 1, 1, 


4_ r9 1 9 9i 
+ [z, l,z,2J 






_Lr9 9 
+ [Z, z, 


1 

i , 


ij 














^'nnn^nnnn9nnn^9nnnn/l^ 

^^U, U, U, U, U, U, U, U, Z, U, U, U, U, Z, U, U, U, U, 


l^ 1 1 


1 
1 


-1- n 1 


9 1 
Z, i 


4- n 


9 19 
Z, i, Z 


4- r9 1 1 1 

+ [z, i, 1, i. 


4- f9 1 1 91 
+ [z, i, 1, ZJ 






4- [9 9 
T[Z, Z, 


9 


11 












<-^64 


(Tinnnnnnnnni 1 nnnnnn^^ 

^U,U,U,U,U,U,U,U,U,U, i, i,U,U,U,U,U,U,4j 


ri 1 1 


1 
1 


J- ri 1 

+ [J-, J- 


9 9 
Z, Z 


4- fl 


9 1 1 

Z, i, i 


4- fl 9 1 9 
+ [J-, Z, i, Z 


4- fl 9 9 11 
+ [J-, Z, Z, iJ 






_i_r9 1 


1 

i , 


11 -1- [9 


1 1 


91 4- 


9 19 

z, J-, z, 


91 1 ro 9 
ZJ + [Z, z, z. 


1 1 1 [9 9 9 91 
iJ + [Z, Z, Z, ZJ 


<-^65 


(Ti nnnnnnnni nni nnnnnA"> 

^U,U,U,U,U,U,U,U,U, i,U,U, i,U,U,U,U,U,4j 


n 1 1 


1 
i 




1 9 


4- n 


19 9 

i, z, z 


4- fl 9 1 9 
[i, Z. i, z 


4- fl 9 9 11 

+ [J-, z, Z, iJ 






4- [9 1 


1 

1 , 


11 -1- f9 
-LJ + [z, 


1 9 

-L, Z, 


11 4- 


9 9 1 

, , , 


91 ^ f9 9 9 
ZJ 1 [z, z, z, 


1 1 ^ f9 9 9 91 




((\ n n fi n 1 fi n fi n 1 n fi n n zt"; 




1 

1. 




9 1 


4- fl 


19 9 
1, Z, Z 


4- fl 9 1 9 
+ [1, Z, i, z 


4- fl 9 9 11 
+ [i, Z, Z, iJ 






4- [9 1 


1 

L . 


11 -1- f9 


1 9 

J-, z, 


91 4- 
ZJ + 


9 9 1 

^, z, i, 


11 4- r9 9 1 

iJ -|- [Z, Z, 1, 


91 4- r9 9 9 91 
ZJ -|- [Z, Z, Z, ZJ 


<-^67 


/"nnnnnnnnnnnnnnnnnnA'i 

U, U, U, U, U, U, U, U, U, U, U, U, U, U, U, U, U, 4J 


n 1 1 


1 
1 


4- n 1 
+ [J-, J- 


1 9 


4- n 


19 1 

1, Z, 1 


fl 9 1 1 

[J-, z, i, 1 


4- fl 9 1 9l 
+ [J-, Z, i, ZJ 






[9 1 

— 1 Z , 1 , 


1 

1 , 


11 4- [9 
J-J + [z, 


1 9 


11 4- 


9 9 1 

z, z, 1, 


1 1 4- [9 9 9 

IJ + [Z, z, z, 


91 




/'nnnnnnnni 1 01 oonnnn^^ 

^U,U,U,U,U,U,U,U, i, i,Z, i,Z,Z,U,U,U,U,oj 


n 1 1 


1 
i 


^ n 1 
+ [1, 1 


9 1 

z, i 


4_ n 
+ [1 


9 1 1 

z, J-, J-, 


4_ [9 1 1 1 

+ [z, 1, i, 1, 


4_ [9 9 9 91 
+ [z, z, Z, /iJ 


(-^69 


(nc\c\(^(\(^c\n^ 91 91 9nnnn'?"\ 
^^u, u, u, u, u, u, u, u, 1, z, 1, z, 1, z, u, u, U, U, oj 


n 1 1 


1 
i 


4- n 1 


1 9 


4- n 


9 1 1 
Z, i, i 


4- f9 1 1 1 

+ [z, i, 1, i. 


4- f9 9 9 91 
+ [z, Z, Z, ZJ 




f'nnnnnnnn9i 1 991 nnnn'?"\ 


n 1 1 
[J-, J-, J- 


1 
± 


4- n 1 
+ [J-, J- 


1 9 
1, Z 


4- fl 


19 1 

1, Z, 1 


4- f9 1 1 1 

+ [z, J-, J-, J-, 


4- f9 9 9 91 
+ [Z, Z, Z, ZJ 




f'nnnnnnnn999i 1 1 nnnn*?"! 


n 1 1 


1 
L 


4- n 1 


1 9 


4- fl 


19 1 

1, Z, 1 


4- fl 9 1 1 

+ [J-, z, i, 1 


4- f9 9 9 91 
+ [z, Z, Z, ZJ 




/'nnnnnnnni 1 1 1 1 1 nnnnQ^ 

^U,U,U,U,U,U,U,U, i, i, i, i, i, i,U,U,U,U,oj 


rill 

[i, i, 1 


1 
i 


^[19 


9 9 

z, z 


4_ r9 
+ [2 


119 

i, i, z 


4_ [9 1 9 1 

+ [z, i, Z, i 


4_ [9 9 1 11 
+ [z, Z, i, iJ 




^U,U,U,U,U,U,U,U,U,U,i,Z,U,U,U,U,U,U,oj 


\^ 1 1 


1 
i 


^ fl 1 


9 9 

z, z 


4_ fl 
+ [1 


9 1 1 

z, i, i 


4_ [9 1 9 1 

+ [z, i, Z, i 


4_ [9 9 1 91 
+ [z, Z, i, ZJ 


O74 


(^U,U,U,U,U,U,U,U,U,U,Z,i,U,U,U,U,U,U,oj 


n 1 1 


1 


^[11 

+ [1, 1 


9 9 

z, z 


4_ n 
+ [1 


9 19 

z, i, z 


1 [9 1 1 1 

+ [z, i, i, 1 


1 [9 9 9 11 

+ [z, /!, Z, iJ 




/'n nnnnnnnni nn9nnnnnQ^ 

(^U,U,U,U,U,U,U,U,U, i,U,U,Z,U,U,U,U,U,oj 


fi 1 1 


1 
1 


4_ n 1 
+ 11, 1 


9 9 

z, z 


4_ n 
+ [1 


9 1 1 

z, i, J-, 


4_ [9 1 1 9 

+ [z, 1, i, z 


4_ [9 9 9 11 
+ [z, z, Z, iJ 


(^76 


^u, u, u, u, u, u, u, u, u, z, u, U, i, u, u, u, U, U, oj 


ri 1 1 


1 
i 


4. ri 1 
+ [1, 1 


9 9' 


4_ fl 
+ [1 


9 9 1' 

z,z, 1 


4_ [9 1 1 1 ' 

+ [z, 1, 1, 1, 


4_ f9 9 1 9l 


C77 


(0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 2, 0, 0, 0, 0, 3) 


[1,1,1 


1 


+ [1,1 


2,1 


+ [1 


2,1,2; 


+ [2,1,1,2; 


+ [2,2,2,1] 


C78 


(0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 1, 0, 0, 0, 0, 3) 


[1,1,1 


1 


+ [1,2 


1,2 


+ [1 


2,2,1; 


+ [2,1,1,1 


+ [2,1,2,2] 


C79 


(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 3) 


[1,1,1 


1 


+ [1,1 


1,2] 


+ [1 


2,2,1] 


+ [2,1,1,1] 


+ [2,1,2,2] 






+ [2,2, 


1, 


2] 












Cso 


(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 3) 


[1,1,1 


1 


+ [1,1 


2,1] 


+ [1 


2,1,2] 


+ [2,1,1,1] 


+ [2,1,2,2] 






+ [2,2, 


2, 


1] 












Csi 


(0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 3) 


[1,1,1 


1 


+ [1,1 


2,2] 


+ [1 


2,1,1] 


+ [2,1,1,1] 


+ [2,2,1,2] 






+ [2,2, 


2, 


1] 














(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3) 


[1,1,1 


1 


+ [1,1 


2,2] 


+ [1 


2,1,1] 


+ [1,2,1,2] 


+ [2,1,1,1] 






+ [2,1,2, 


l] + [2,2,2. 


2] 
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TABLE XI. Representative elements for the sets of equivalence classes for n = 4, Z) = (2, 2, 2, 2) 
induced by the permutation symmetry of the spaces in {Vi, . . . , V4}. A representative element v 
is given hy v = Avi, where A: V V is sl certain linear operator and G F is a fixed vector. 
(Without loss of generality and for comparison with other tables, we choose vi = [1, 1, 1, 1].) The 
operator a„j is defined by aj[. . . , 1, . . . ] = [. . . , 2, . . . ] and aj[. . . , 2, . . . ] = [...,!,...], where only 
the ith index changes. To obtain all classes in each group, all possible choices of the indices 
{i,j,k,l} = {1,2,3,4} should be considered. This results in 27 fundamental sets of 83 classes. 



Co 

Ci 

{Cs, Cg, Cio} 

{Cii, . . . , C14} 

{Ci5, . . . , Cis} 

C'lg 

{C20, • • • , C25} 
C26 

{C27, C28, C29} 
{C30, C31, C32} 
C33 

{C34, . . . , C37} 
{C38, C39, C40} 
{C41, . . . , C46} 
{C47, C48, C49} 
{C50, . . . , C53} 
{C54, . . . , C59} 
C'eo 

{Cei, C62, Ces} 
{C64, Ces, Cee} 
Cer 

{Ces, • • • , C71} 
C72 

{C73, . . . , C78} 
{C79, Cgo, Csi} 

C^2 





1 

1 + 

(n 
1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 - 
1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 - 
1 + 
1 + 
1 + 
1 + 
1 



- aiaj){l + Qkm) 
ai{aj + Uk) 

ai{aj +ak + ai) 
ai{aj + Qkai) 
QiajUkai 

(ttj + aj){ak + ai) 

[ai + aj + ttfe + ai) + ajakai + atakai + aiUjai + aiajak + Ojaj + a^a; + aiGjUkai 
c{aiaj + akai) - (1 + c){aiak + ajai) + aiajakai, cG F, {-2, -1,0, 1} 
Gi + ajak + aiUkai 
aittj + atai 
{ai + aj){l + akai) 
aiUj + (flj + aj)akai 
ai{aj + ttk) + ajttkai 
ataj + akai + a^ak 
{ai + aj){ak + a;) + aiajakai 
a 



+ aj + akai + aiak + aiajai 
+ Gj + aktti + aiai + ajUk + ajai + aiGkai + atajak + aiajakai 
- aj +ak + ai + Uiaj + aiak + ajai + aiajakai 

+ aj + ttfe + a^ajakai 
{aj + flfe + ai) + ajakai 
(1 + akai) + aj{ak + ai) 
+ Oj + akai + aiaj{ak + at) 



+ (I, + (I j + <ll:<ll + <<i<<k + "j"' + "i" 



V. COMPARISON WITH CLASSICAL INVARIANTS 



The central distinction between the classification presented in this work and classifications 
found in the literature is in the type of invariants used. We rely on discrete invariants, while 
most other methods use continuous invariants. Broadly speaking, the relation between the 
two types of invariants is such that different values of the discrete invariants correspond to 
certain continuous invariants being zero or nonzero. We do not attempt here the complete 
comparison between classifications based on the two types of invariants and present results 
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only for the methods reviewed in [19] for 3 qubits and developed in [8] for 4 qubits. 
A. n — 3 

For an arbitrary state of 3 qubits 

2 2 2 

jl=l j2 = l j3 = l 

the classical invariants (see e.g. [19]) are 

= t'1,1,11'1,2,2 - t'l,l,2t'l,2,l + t'2,l,lt'2,2,2 - t'2,l,2t'2,2,l , 
^2 = Vx,\,\V2,\,2 - Vx,\,2V2,\,\ + ^'l,2,l^'2,2,2 " ^'l,2,2t'2,2,l , 

^3 = ^'l,l,l^'2,2,l - ^'l,2,l^'2,l,l + ^'l,l,2^'2,2,2 " ^'l,2,2t'2,l,2, 

7 2 2,2 2,2 2,2 2 

Al4 — ^1^1^1^2,2,2 "T '^1,1,2'^2,2,1 + '^l,2,l'^2,l,2 + '^2,l,l'^l,2,2 

- 2(f i,i,it'i,i,2^^2,2,1^^2,2,2 + ^^1,1,1^^1,1,2^'2,2,1^'2,2,2 + ^'l,l,l^'l,l,2^^2,2,l^'2,2,2 

+ ■1^1,1,11^1,1,21^2,2,1^2,2,2 + ■J^l,l,l,l'J^l,l,2l^2,2,l1^2,2,2 + 1^1,l,l1^1,l,2l^2,2,l1^2,2,2) 

+ 4(t'i,i,it'i,2,2'y2,l,2'y2,2,l + ^'l,l,2^'l,2,l^'2,l,l^'2,2,2)- 

Table XII lists zero and nonzero values of /ii, . . . , /i4 for the classes Cq, . . . , Ce- Ignoring 
the trivial difference between Cq and Ci, we conclude that zero values of the continuous 
invariants hi, . . . ,h4 distinguish all the classes for 3 qubits found using the discrete invariants. 
Similar comparisons for n = 3, D = {di, d2, ds) can be easily done. 



TABLE XII. The continuous invariants /ii, . . . , /i4 for the sets of equivalence classes based on the 
discrete invariants for n = 3, Z) = (2, 2, 2). 





hi 


h2 


h3 


/14 


Co 














Ci 






















7^0 





Cs 





7^0 








Ci 


7^0 











C5 


7^0 


7^0 


7^0 





Ce 











7^0 



B. n = 4 

For an arbitrary state of 4 qubits 

2 2 

Jl=l j4=l 
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the Hilbert series lead [8] to the classical invariants 



ho 



/l5 
^5,1,1 
^5,1,2 
^5,1,3 
^5,2,1 
^5,2,2 



^5,2,3 = 
^5,3,1 = 
^5,3,2 = 
^5,3,3 — 

he = 

^6,1,1 = 
^6,1,2 = 
^6,1,3 = 

^6,2,1 = 



6,2,2 



6,2,3 



h 

^6,3,1 
^6,3,2 
^6,3,3 

hr 

^7,1,1 



'^l,l,l,l'i^2,2,2,2 ~ 'i^l,l,l,2'i^2,2,2,l ~ 'i'l,l,2,l'i'2,2,l,2 + 'i'l,l,2,2'i'2,2,l,l 
■1^1,2, l,l1'2,l,2,2 + 1'l,2,l,2l'2,l,2,l + 1'l,2,2,l1'2,l,l,2 — 1'l,2,2,2l'2,l,l,l , 

^1,1,1,1 ^'1,2,1,1 ^'2,1,1,1 ^^2,2,1,1 

"^1,1,1,2 '^1,2,1,2 '^^2,1, 1,2 '^2,2,1,2 

t'1,1,2,1 ^"1,2,2,1 ^^2,1,2,1 ^"2,2,2,1 

■J^l, 1,2,2 '^^1,2,2,2 '1^2,1,2,2 '^'2,2,2,2 

^^l,l,l,l ^^2,1,1,1 -^1,1,2,1 ^^2,1,2,1 

^'1,1,1,2 '^2,1,1,2 '^1,1,2,2 '^2,1,2,2 

""1,2,1,1 ^2,2,1,1 ^"1,2,2,1 ^"2,2,2,1 

^"1,2,1,2 ""2,2,1,2 ""1,2,2,2 ""2,2,2,2 

^^1,1,1,1 ^^1,1,1,2 ^^2,l,l,l ^^2,1,1,2 
""1,1,2,1 ""1,1,2,2 ""2,1,2,1 ""2,1,2,2 
""1,2,1,1 ""1,2,1,2 ""2,2,1,1 ""2,2,1,2 
""1,2,2,1 ""1,2,2,2 ""2,2,2,1 ""2,2,2,2 
= dct(/?.5 jj)i<j<3.i<j<3, 
= -?^1,1,1,2?^1, 1,2,1 + ^^1, 1,1, 1^1, 1,2,2, 

= '"l,l,2,2'"l,2,l,l — ^l,l,2,l""l,2,l,2 — ""l,l,l,2""l,2,2,l + ""l,l,l,l""l,2,2,2, 
= ~""l,2,l,2'"l,2,2,l + ""l,2,l,l""l,2,2,2; 

= f l,l,2,2t'2,l,l,l - t'l,l,2,lt'2,l,l,2 - t'l,l,l,2t'2,l,2,l + t'l,l,l,lt'2,l,2,2, 
= ""l,2,2,2""2,l,l,l ~ ""l,2,2,l""2,l,l,2 ~ ""l,2,l,2""2,l,2,l + ""l,2,l,l""2,l,2,2 
+ ""l,l,2,2""2,2,l,l — ""l,l,2,l""2,2,l,2 — ""l,l,l,2""2,2,2,l + ""l,l,l,l""2,2,2,2, 
— ""l,2,2,2""2,2,l,l ~ ""l,2,2,l""2,2,l,2 ~ ""l,2,l,2""2,2,2,l + ""l,2,l,l""2,2,2,2, 
= — ""2,1,1,2""2,1,2,1 + ""2,1,1,1""2,1,2,2, 

""2,1,2,2""2,2,1,1 ~ ""2,l,2,l""2,2,l,2 ~ ""2,1,1,2""2,2,2,1 + ""2,l,l,l""2,2,2,2, 
~""2,2,1,2""2,2,2,1 + ""2,2,1, l""2,2,2,2, 
det(/l6,i,j)l<i<3;l<j<3, 
-^^1,1,1,2^1,2,1,1 + ^^1,1,1,1^1,2,1,2, 

-""1,1,2,2""1,2,1,1 + ""1,1,2,1""1,2,1,2 - ""1,1,1,2^^1,2,2,1 + ""1,1,1,1^^1,2,2,2, 
— ""l,l,2,2""l,2,2,l + ""l,l,2,l""l,2,2,2, 

^"1,2,1,2^2,1,1,1 - ""l,2,l,lt'2,l,l,2 - ""l,l,l,2t'2,2,l,l + ""l,l,l,l^'2,2,l,2, 

""1,2, 2,2""2, 1,1,1 — "?-'l,2,2,l""2,l,l,2 + "?-'l,2,l,2"'-'2,l,2,l "?-'l,2,l,l"'-'2,l,2,2 
""l,l,2,2""2,2,l,l + '"l,l,2,l""2,2,l,2 — ""1,1,1,2^2,2,2,1 + '"l,l,l,l'"2,2,2,2, 
""1,2, 2,2""2, 1,2,1 ~ ""l,2,2,l""2,l,2,2 — ""l,l,2,2""2,2,2,l + ""l,l,2,l""2,2,2,2, 
~""2,1,1,2""2,2,1,1 + ""2,1,1,1""2,2,1,2, 

— ""2,1,2, 2'"2,2, 1,1 + ""2,1,2,1""2,2,1,2 — ""2,1,1,2""2,2,2,1 + ""2,l,l,l""2,2,2,2 , 
— ""2,1,2, 2""2,2, 2,1 + ""2,1, 2,1""2, 2,2,2, 
det (/t7,j,j ) l<i<3;l<j<3 , 
-""1,1,2,1^1,2,1,1 +'"1,1,1,1'"1,2,2,1, 
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^7,1,2 = -^^1,1,2,2^1,2,1,1 - ^^1,1,2,1^^1,2,1,2 + ^1,1,1,2^1,2,2,1 + ^1,1,1,1^1,2,2,2, 

^7,1,3 — — 1'l,l,2,2'J^l,2,l,2 + 'J^1,1,1,2'J^1,2,2,2, 

^7,2,1 = ■J^1,2,2,1'J^2,1,1,1 - ■J^1,2,1,1'J^2,1,2,1 " 'yi,l,2,l'J^2,2,l,l + ■J^l,l,l,l'J^2,2,2,l , 

^7,2,2 — 'l'l,2,2,2'l'2,l,l,l + 'f^l,2,2,l'f^2,l,l,2 ~ 'l'l,2,l,2'l'2,l,2,l ~ 'l'l,2,l,l'l'2,l,2,2 

— 'yi,l,2,2'^2, 2,1,1 — '^l,l,2,l'y2,2,l,2 + 1^1,l,l,2'i^2,2,2,l + '^l,l,l,l'i'2,2,2,2, 

^7,2,3 — 'i^l,2,2,2'y2,l,l,2 — 'yi,2,l,2'y2,l,2,2 " t^l, 1,2,2^2,2,1,2 + 1^1,1, l,2'y2,2,2,2, 

^7,3,1 = ~'y2,l,2,l'i^2,2,l,l + '^^2, 1,1, l'i^2, 2,2,1; 

^7,3,2 = — 'i^2,l,2,2'i^2,2,l,l — 'i^2,l,2,l''-'2,2,l,2 + 1^2,1, 1,21^2,2,2,1 + 1^2,1, 1,11^2,2,2,2, 

^7,3,3 = ~1^2,l,2,2l^2,2,l,2 + 1^2,1,1,21^2,2,2,2- 

The invariants satisfy the relations 

^2 + ^3 + ^4 = 0, 
^1^2 - /l6 + ^7 = 0, 

hihs — + = Q, 
hih^ — + he — 0, 

which imply that only two invariants among ^2, ^3, ^4 are independent and only one invariant 
among htj,hQ, h-j is independent. Nevertheless, in the results below we use all the invariants 
/ii , . . . , /i7 for the sake of symmetry. 

Tables XIII, XIV, XV hst zero and nonzero values of /ii, . . . , /ly for the classes Co, ... , €^2- 
Tabic XVI shows that with only 21 possibilities for their independent values, the zero values 
of invariants /ii, . . . , /17 cannot distinguish all the classes found using the discrete invariants. 
Of course, since hi,...,h7 is a complete list of invariants for 4 qubits, all their possible 
values completely characterize the entangled states (and with a greater refinement than our 
invariants, as explained earlier). The difficulty of using the continuous invariants is of course 
in finding all their allowed values. 

Note that some of our classes split with respect to the values of the classical invariants, 
but it should be remembered that there are relations that could allow individual classical 
invariants to be less constrained than their irreducible set. 

We also note that the families of entangled 4 qubits found in [5] are related to our 
entanglement classes as follows: 

Gabcd G Cs2, Labc2 ^ Cs2, -^'02 62 ^ {^"73, ■ ■ ■ , Cyg}, Lab^ € Cqt, 

La4 £ {C73, ■ ■ ■ , Cys}, -^0203^1 e {Ces, ■ ■ ■ , Cji}, -^05^3 £ {C34, ■ ■ ■ , Csr}, 

^Orsn ^ {^50, ■ ■ ■ , ^53}, Log^j-Ogg^j- e {Ci5, . . . , Cis}. 

With the correction suggested in [18], we find the same result as above except now Lab^ G Cs2- 
This leaves many classes found here unaccounted for in the method of [5] . 

VI. CONCLUSIONS 

Mathematical structure of entangled states gives rise to new entanglement invariants, 
which lead to a new method of entanglement classification. The invariants describe alge- 
braic properties of linear maps associated with the states. For finite-dimensional spaces. 
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TABLE XIII. The continuous invariants hi,. . . ,hf for the sets of equivalence classes based on the 
discrete invariants for n = 4, D = (2, 2, 2, 2). 







"•2 


"•3 


1, 

/14 


u 

"5 




1, 

"7 




u 


n 

u 




u 






u 




n 
u 


n 
U 


n 
u 


n 
u 


n 
u 


u 


fl 




u 


u 


u 


u 


u 


u 


u 


<^3 


u 


n 
u 


u 


u 


u 


u 


fl 


O4 


u 


u 


u 


u 


u 


u 


n 
u 




u 


n 
u 


u 


u 


u 


u 


n 
u 




u 


u 


u 


u 


u 


u 


n 
u 




n 
u 


n 
u 


u 


u 


u 


u 


fl 


ri 






-J- n 
7^ u 


n 

u 


U 


u 


-J- n 




7^ u 




n 
U 


-J- n 
7^ U 


n 
U 




U 




7^ u 


U 




7^ u 




n 
U 


U 


Oil 


U 


U 


U 


U 


U 


fl 
U 


fl 
U 


O12 


u 


u 


u 


u 


u 


fl 
U 


fl 


L-ia 


u 


n 


u 


u 


u 


fl 

u 


fl 
u 


Oi4 


u 


n 
u 


u 


u 


u 


fl 
u 


fl 
u 


<-^15 


u 


U 


u 


u 


u 


fl 

u 


fl 
u 


<-^16 


u 


u 


u 


u 


u 


fl 
u 


fl 
u 


{_/17 


u 


u 


u 


u 


u 


fl 
u 


fl 
u 


<-l8 


u 


u 


u 


u 


u 


fl 
u 


fl 
u 


<-^19 


u 


u 


u 


u 


u 


fl 
u 


fl 
u 




u 


n 
u 


u 


u 


u 


fl 
u 


n 
u 


O2I 


u 


u 


u 


u 


u 


fl 
u 


fl 




u 


n 

u 


u 


u 


u 


fl 

u 


fl 


C23 


u 


n 

u 


n 
u 


u 


u 


fl 
u 


fl 


C24 


n 
u 


n 


n 


n 
w 


n 


fl 


n 
u 


/-^ 
<-^25 


n 


n 


n 


n 


n 


fl 


n 
u 


/-^ 

'-^26 




n 

u 


n 
u 


n 
u 


fl 
u 


fl 
u 


n 
u 


O27 




7^ ^ 




U 


fl 
U 


fl 

u 


fl 

7^ u 


<-^28 


7^ u 


7^ u 


u 


7^ u 


fl 

u 


-4- fl 


fl 
u 


(-^29 




n 
u 




7^ u 




fl 

u 


fl 

u 


<-30 


7^ u 






u 




-1- n 


-A fl 


<-^31 




7^ u 


U 




fl 

7^ u 


-4- n 


fl 


<-^32 




u 


7^ u 




-A fl 


-4- n 


fl 

7^ u 




5^ 


^ 


5^ 


^ 


9^ 


5^ 


9^ 


C34 























C35 























C36 























C37 























C38 


7^0 




















C39 


7^0 




















C40 


7^0 
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TABLE XIV. The continuous invariants hi,. . . ,hf for the sets of equivalence classes based on the 
discrete invariants for n = 4, D = (2, 2, 2, 2). 







"2 


113 


A14 




n% 


rir 


041 








n 
u 


u 


u 


-J- n 


(-^42 




-J- n 


7^ U 


n 
U 


n 
U 


u 


n 

T U 


(-^43 




^ n 
7^ u 


U 


7^ ^ 


n 
u 


7^ L) 


U 


O44 




^ n 


n 

u 


-J- n 
7^ U 


n 

u 


7^ U 


U 


L'45 




n 
U 


7^ U 


-4- n 


-4- n 


U 


U 


C4(i 




n 
U 








U 


n 
U 


O47 


U 


7^ 


7^ 


U 


u 


U 


U 


r" 


7^ u 




7^ 


n 
(J 


(J 


(J 


7^ U 


(-^48 


u 


7^ U 


n 

U 




u 


U 


u 


•^48 




7^ U 


u 


7^ u 


u 


7^ U 


u 


<-^49 
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TABLE XV. The continuous invariants /ii , . . . , /17 for the sets of equivalence classes based on the 
discrete invariants for n = 4, D = (2, 2, 2, 2). 
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each invariant takes a value from a finite set of integers, and the resulting set of entan- 
glement classes is finite. The relation to the standard continuous invariants is such that 
different values of the discrete invariants correspond to certain continuous invariants being 
zero or nonzero. We believe that our classification is the most refined restricted classification 
possible. Although this result is formulated as a theorem in the text, its proof is not a usual 
mathematical proof, but rather a proof by exhaustion of possibilities. 

The new method works for an arbitrary finite number of spaces of finite dimensions. As 
its application, we obtained entanglement classifications for a wide selection of individual 
cases of three subsystems and the case of four qubits. 

For three subsystems, in addition to finding classifications for individual values oi D = 
{di, d2, d), it is rather easy to obtain results for infinite sequences of values of d. An inter- 
esting general feature of these results (which is easy to prove) is that increasing d beyond 
did2 does not introduce any new entanglement classes. As examples, we have found such 
classifications for the values ((ii^fia) G {(2, 2), (2, 3), (2, 4), (2, 5), (2, 6), (3, 3), (3, 4)} and ar- 
bitrary d. Only one of these sequences, D = (2, 2, d), had been conjectured in the literature, 
for which our method gives the same number of classes as the classifications in [3], [4], [7], 
[9] and the conjectured classification in [7], [9]. 

Entanglement classes and representative elements could be easily generated for other 
infinite sequences. The classification problem for the general case of three subsystems, 
however, is challenging and currently under study. Note that the entanglement of a set of 
three large spin subsystems is in some practical sense complementary to a system of many 
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TABLE XVI. The continuous invariants hi,. . . ,h'r for the sets of equivalence classes based on the 
discrete invariants for n = 4, D = (2, 2, 2, 2). 
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low spin (e.g., many qubit) subsystems. Both have potential for the construction of practical 
devices. 

The classification of entanglement of four qubits has been considered by several groups 
of authors [5, 13-18, 20]. All previous works found 9 or fewer fundamental sets of classes 
after permutations have been removed. In our work wc found 27 fundamental sets of classes. 
Our refined classification could be useful to experimenters who consider detailed properties 
of four qubit systems. For example, Barreiro et al. [23] find a rich dynamics when they 
arrange four Ca^ ions as qubits and study entanglement via decoherence and dissipation. 
Sec also [24] for earlier 4 qubit work. 

To deepen our knowledge about other quantum systems, their entanglement should be 
thoroughly studied as well. Our method provides a simple, general, practical approach to 
such studies. 

Our new invariants are topological since they are the dimensions of linear spaces. Al- 
though the invariants are rather simple from the point of view of topology, they may have 
a different interpretation when viewed from other perspective. For an example of a possibly 
related interpretation, see [25]. Finally, it is also worth pointing out that while we find pure 
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representative states for each class, it is straightforward to combine them into mixed states 
via a density matrix approach. 
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